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The hydrodynamic slippage at a solid-liquid interface is currently at the center of our un- 
i-^j , derstanding of fluid mechanics. For hundreds of years this science has relied upon no- 

slip boundary conditions at the solid-liquid interface that has been applied successfully to 
model many macroscopic experiments, and the state of this interface has played a minor 
O . role in determining the flow. However, the problem is not that simple and has been revis- 

ited recently. Due to the change in the properties of the interface, such as wettability and 
| roughness, this classical boundary condition could be violated, leading to a hydrodynamic 

slip. In this chapter, we review recent advances in the understanding and expectations for 
' the hydrodynamic boundary conditions in different situations, by focussing mostly on key 

, papers from past decade. We highlight mostly the impact of hydrophobicity, roughness, 

\Q • and especially their combination on the flow properties. In particular, we show that hy- 

l/"") ' drophobic slippage can be dramatically affected by the presence of roughness, by inducing 

^v^j , novel hydrodynamic phenomena, such as giant interracial slip, superfluidity, mixing, and 

low hydrodynamic drag. Promising directions for further research are also discussed. 

1.1. Introduction 

^ ' Fluid mechanics is one of the oldest and useful of the 'exact' sciences. For hun- 

■ dreds of years it has relied upon the no-slip boundary condition at a solid-liquid 

interface, that was applied successfully to model many macroscopic experiments.^ 
However, the problem is not that simple and has been revisited during recent 
years. One reason for such a strong interest in 'old' problem is purely fundamen- 
tal. The no-slip boundary condition is an assumption that cannot be derived from 
first principles even for a molecularly smooth hydrophilic [the contact angle (fixed 
by the chemical nature of a solid) lies between 0° and 90°] surface. Therefore, the 
success of no-slip postulate may not always reflect its accuracy but in fact rather 
the insensitivity of the experiment. Another reason for current interest to flow 
boundary conditions lies in the potential applications in many areas of engineer- 
ing and applied science, which deal with small size systems, including micro- and 
nanofluidics, 2 flow in porous media, friction and lubrication, and biological flu- 
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ids. The driving and mixing of liquids when the channel size decreases represent a 
very difficult problemP There is therefore a big hope to cause changes in hydrody- 
namic behavior by an impact of interfacial phenomena on the flow. For example, 
even ideal solids, which are both flat and chemically homogeneous, can have a 
contact angle, which exceeds 90° (the hydrophobic case). This can modify hydro- 
dynamic boundary conditions, as it has been shown yet in early work. 4 Besides 
that, solids are not ideal, yet rough. This can further change, and quite dramati- 
cally, boundary conditions. It is of course interesting and useful to show how the 
defects or pores of the solids modify them. But today, the question has slightly 
shifted. Thanks to techniques coming from microelectronics, we are able to elab- 
orate substrates whose surfaces are patterned (often at the micro- and nanometer 
scale) in a very well controlled way (see Fig. Il.l.jl , which provides properties (e.g. 
optical or electrical) that the solid did not have when flat or slightly disordered. 
A texture affects the wettability and boundary conditions on a substrate, and can 
induce unique properties that the material could not have without these micro- 
and nanostructures. In particular, in case of super-hydrophobic solids, which are 
generated by a combination of surface chemistry and patterns, roughness can dra- 
matically lower the ability of drops to stick, by leading to the remarkable mobility 
of liquids. At the macroscopic scale this renders them 'self-cleaning' and causes 
droplets to roll (rather than slide) under gravity and rebound (rather than spread) 
upon impact instead of spreading.^ At the smaller scale, reduced wall friction and 
a superlubricating potential are almost likely associated with the breakdown of the 
no-slip hypothesis. This suggests that super-hydrophobic surfaces should strongly 
affect the transport of fluids. This area of research is relatively new, but rapidly de- 
veloping, and it already attracted scientists from physics, chemistry, mathematics, 
and engineering. 




Figure 1.1. Examples of anisotropic (grooves) and isotropic (pillars, chess-board) textures. 

In this chapter we review recent advances in the understanding and expec- 
tations for the fluid-solid boundary conditions in different situations, where hy- 
drophobicity and roughness impact the flow properties, by focussing mostly on 
key papers from past decade. Throughout the chapter, we emphasize open ques- 
tions. We start with a very brief history of the no-slip boundary conditions. Then, 
after introducing the terminology and models, and describing new developments 
and instruments, which give the possibility of investigating fluid behavior at the 
micro- and nanoscale, in the following section we present results obtained for 
smooth surfaces, by highlighting the role of wettability. Then follows the results 
for a rough hydrophilic and, especially, hydrophobic surfaces. In the latter case we 
show how can roughness enhance hydrodynamic slip and thus the efficiency of 
interfacial transport phenomena. We close with the suggestions to guide the next 
round of experimental and theoretical studies in this expanding area of research. 
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1.2. Origins and history 

The nature of boundary conditions in hydrodynamics was widely debated in 19th 
century and many of the great names in fluid dynamics have expressed an opinion 
on the subject, as discussed in a recent reviewP The linear slip boundary condi- 
tion was introduced by Navierj^ and this remains a standard characteristic of slip 
used today. Helmholtz and von Piotrowski were probably the first to report some 
evidence of slippage at the solid-liquid interface. 8 We refer the reader to a compre- 
hensive review 9 for detailed accounts of early experimental work. The significance 
of hydrophobicity for the slip phenomenon does not seem to be recognized at this 
early stage. 

The matter was revisited at the beginning of 20th century as reviewed in SI Sev- 
eral research groups reported faster flow after they had treated the capillary tube 
with organic compounds, and interpreted results in favour of slippage of water. 
There have been however some conclusions about complications from surface ten- 
sion effects in these early experiments. Although not recognized by the authors, 
their work tried to make a connection between the amount of slip and the hy- 
drophobization of the solid. That is why, despite the fact that the early results still 
did not provide a strong support to the hypothesis of a hydrophobic slippage, we 
consider these papers to be an important contribution to the subject. The first re- 
liable results for water were probably obtained at the second part of 20th century. 
The measurements by Schnell 10 and Churaev et alPI in thin hydrophobic capillar- 
ies unambiguously suggested a co ncept o f hydrophobic slippage, which allowed 
to formulate first theoretical models!- 4 ! 12 ! 13 ! However, despite increasing body of ev- 
idence in favor of hydrophobic slippage, experimental methods and results have 
been rather limited, the amplitude of slip and its dependence on various parame- 
ters (such as the contact angle, shear-rate, etc) still remained an open question. 

In 21st century, or during the last decade, the field of hydrodynamic boundary 
conditions is rapidly advancing, motivated by potential applications in microflu- 
idics, as well as fundamental scientific issues in colloid physics and fluid mechan- 
ics. Both hydrophobic and rough interfaces have received enormous attention. 
Note that some of the first data completely escaped from the theoretical frame- 
work, with qualitative (shear rate threshold to slippage; shear rate dependent slip 
length) and quantitative (micrometric versus nanometric slip length) discrepan- 
cies,^ which have caused hot debates. Fortunately, an avalanche development of 
the field, including experimental methods, allowed one, and relatively fast, clar- 
ify the situation and highlight reasons for existing controversies. Thus, it is pos- 
sible to conclude that boundary conditions at the smooth hydrophobic surfaces 
are currently quite well understood. We also suggest that a rough interface has 
recently been relatively well understood, although there are still some open ques- 
tions. Over the past few years hot and promising directions of research seem to be 
the strategy of a combination of hydrophobic slip with surface roughness and/ or 
interfacially driven flow, such as electro-osmosis. These new systems of interest 
could potentially lead to a generation of extremely fast flows in microfluidic de- 
vices and stimulated new fundamental studiesfSI 
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1.3. Terminology and models 

We will refer to as a slip any situation where the value of the tangential compo- 
nent of velocity appears to be different from that of the solid surface. The simplest 
possible relation assumes that the tangential force per unit area exerted on the 
solid surface is proportional to the slip velocity. Combining this with the consti- 
tutive equation for the bulk Newtonian fluid one gets the so-called (scalar) Navier 
boundary condition^ 

u s = b^ r (1.1) 
dz 

where n s is the (tangential) slip velocity at the wall, du/dz the local shear rate, and 
b the slip length. This slip length represents a distance inside the solid to which the 
velocity has to be extrapolated to reach zero. The standard no-slip boundary con- 
dition corresponds to b = 0, and the shear-free boundary condition corresponds to 
b — > oc 12 (see Fig. ll.2.t . In the most common situation b is finite (a partial slip) and 
associated with the positive slip velocity. It can, however, be negative, although in 
this case it would not have a lone-ranee effect on the flowP 
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Figure 1.2. Interpretation of the Navier slip length b. 

Obviously the control of slip lengths is of major importance for flow at interface 
and in confined geometry. It would be useful to distinguish between three different 
situations for a boundary slip since the dynamics of fluids at the interface introduce 
various length scales. 

• Molecular (or intrinsic) slip, which allows liquid molecules to slip directly 
over solid surface (Fig. I1.3.h ). Such a situation is not of main concern here since 
molecular slip cannot lead to a large jJMand its calculations requires a molecular 
consideration of the interface region. 

An important point to note that although Eq.(Tj} is the most commonly used 
boundary condition for hydrodynamic slip, it is not widely appreciated thaP also 
postulated the more general relation, Au = Mt where T is the local shear stress 
(normal traction) and M is a constant interfacial mobility (velocity per surface 
stress). For a Newtonian fluid, t = rjdu/dz, this reduces to il.lt with b = Mt], 
where t] is the viscosity. Molecular dynamics (MD) simulations have shown that 
the equation with constant M is more robust than jl.lt with constant b, since the 
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Figure 1.3. Schematic representation of the definition of intrinsic (a), apparent (b), and 
effective (c) slip lengths. Reprinted with permission from!!4l Copyright (2011) by the IOP 
Publishing Ltd. 



fluctuating slip velocity correlates better with the shear stress (normal forces) than 
with velocity gradients very close to the surfaced 

A possible starting point for molecular modeling is a scalar Einstein relation, 
M = bkT I A, relating the M to the 'interfacial diffusivity' per unit area A, by anal- 
ogy with the theory of Brownian motion. This can be recast as a type of a Green- 
Kubo formula, and allows one to prove that the slip length does not depend on the 
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bulk viscosity, rj (since it depends on //D)P§l Rough estimates of the molecular slip 
suggest b oc fcg TrjD/ e 2 , where e reflects a typical energy of interaction between 
liquid and solid.™ By relating e and a contact angle, cos 9 one can then pro- 
pose a scaling relation b oc (1 + cos 6)~ 2 , which is consistent with MD results P3ED 
In other words, the static contact angle is the crucial parameter controlling water 
slippage over smooth hydrophobic surfaces. We remark and stress, however, that 
for realistic contact angles only b below 10 nm has been predicted. Therefore, it is 
impossible to benefit of such a slip in a larger scale applications. 

• The intrinsic boundary condition maybe rather different from what is probed 
in flow experiment at larger length scale. It has been proposed^ to describe the 
interfacial region as a lubricating 'gas film' of thickness e of viscosity rjg different 
from its bulk value rj. A straightforward calculations give apparent slip (Fig. ll.3.b ) 

' (1.2) 



This represents the so-called 'gas cushion model' of hydrophobic slippa gejUwhich 
got a clear microscopic foundation in terms of a prewetting transition!^ Being a 
schematic representation of a depletion close to a wallPthis model provides a use- 
ful insight into the sensitivity of the interfacial transport to the structure of inter- 
face. Similarly, electrokinetic flow displays apparent slip. Note that recent molec- 
ular slip studies™ also suggested a kind of the 'vapor cushion model', where 
b oc e 4 , but apparently in a one-component system the value of e is too small to 
describe most of experimental data, suggesting that a two-component systerrP^QI 
is required. Whether of not the addition of a hydrophobic solute will lead to a non- 
linear dependence of b on e remains an open question and has to be investigated. 

A modification of this scenario would be a nanobubble coated surface, i.e. a 
heterogeneous two-phase depletion layerP2E2 However, as it has been discussed 
yet irp' nanobubbles will reduce the viscosity of the surface layer only if surface 
tension is small enough (which means that the bubbles deform). If the shape of 
the nanobubbles will remain a spherical lens, t]g should exceed the bulk viscosity. 
There remains still an open question connected with the exact expression for the 
effective viscosity of the layer of surface nanobubbles. 

• Another situation is that of effective slip, b e ff , which refers to a situation where 
slippage at a complex heterogeneous surface is evaluated by averaging of a flow 
over the length scale of the experimental configuration (e.g. a channel etc)P^^Ql 
In other words, rather than trying to solve equations of motion at the scale of the 
individual corrugation or pattern, it is appropriate to consider the 'macroscale' 
fluid motion (on the scale larger than the pattern characteristic length or the thick- 
ness of the channel) by using effective boundary conditions that can be applied at 
the imaginary smooth surface. Such an effective condition mimics the actual one 
along the true heterogeneous surface. It fully characterizes the flow at the real sur- 
face and can be used to solve complex hydrodynamic problems without tedious 
calculations. Such an approach is supported by a statistical diffusion arguments 
(being treated as an example of commonly used Onsager-Casimir relations for 
non-equilibrium linear response)/- theory of heterogeneous porous materials,^ 
and has been justified for the case of Stokes flow over a broad class of surfacesPS 
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For anisotropic textures b e ff depends on the flow direction and is generally a ten- 
sor.^ Effective slip also depends on the interplay between typical length scales of 
the system as we will see below. Well-known examples of such a heterogeneous 
system include composite superhydrophobic (Cassie) surfaces, where a gas layer 
stabilized with a rough wall texture (Fig. I1.3.b ). For these surfaces effective slip 
lengths are often very large compared with the value on flat solids, similarly to 
what has been observed for wetting, where the contact angle can be dramatically 
enhanced when surface is rough and heterogeneous.^ 



1.4. Experimental methods 

The experimental challenge generated a considerable progress in experimental 
tools for investigating flow boundary conditions, using the most recent devel- 
opments in optics and scanning probe techniques. A large variability still exists 
in the results of slip experiments so it is important first to consider the different 
experimental methods used to measure slip. Two broad classes of experimental 
approaches have been used so far: indirect and direct (local) methods. 

• High-speed force measurements can be performed with the SFA (surface 
forces apparatus jEKH or AFM (atomic force microscope) PS In particular, in the 
drainage 

methocPESlthe 

end of the spring away from the attached sphere is driven 
toward the (fixed) plane with a constant driving speed (as shown in Fig. II. 4.1 . The 
sphere itself, however, does not move at a constant speed, so that the spring is de- 
flected as a result of both the surface force (which should be measured separately) 
and the hydrodynamic forces. The solution of the (differential) equation of motion 
allows to deduce a drag force, with the subsequent comparison with a theory of a 
film drainage 

1121361 w here the drag force is 

F= _6^| r=fTr (13) 

where — dh/dt is the relative velocity of the surfaces, Fj is the Taylor solution for 
no-slip surfaces (which, however, never appeared in any of G. I. Taylors publica- 
tions as discussed in®]), and /* represents a correction function for slippage, which 
depends on b/h. For example, for a very convenient experimental configuration 
of the interaction of a hydrophilic surface with a hydrophobic one, which allows 
us to avoid the formation of a gas bridge, 37 the correction for slip takes the fornP 



f =T T = l{ 1+ 2b 



ib V h 



(1.4) 



Note that when h <C ib, f* — > 1 /4, so that the force is still inversely proportional 
to the gap, but becomes four times smaller than predicted for hydrophilic surfaces. 
For two similar hydrophobic surfaces, /* — > (h/3b) \n(6b/h) when h -C 6b, so that 
a hydrodynamic drag is only logarithmically dependent on the separation.^ 

Note that AFM force balance incorporates both (concentrated) force on the 
sphere and the drag on the cantilever as shown in Fig. 11.4.1 The drag on a can- 
tilever is neither small nor neglig ible,™ and its ignorance might cause wrong 
experimental conclusions. 
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This approach, being extremely accurate at the nanoscale, does not provide 
visualization of the flow profile, so that these measurements are often identified as 
indirect. 



AFM laser 




Figure 1.4. Schematic of the dynamic AFM force experiment. Adapted frorrP^ 

• Direct approaches to flow profiling, or velocimetry, take advantage of var- 
ious optics to monitor tracer particles. These methods include TIR-FRAP (to- 
tal internal reflection - fluorescence recovery after photo-bleaching)S21 ^-prvlilUS 
(particle image velocimetry), TIRV (total internal reflection velocimetry),* 3 EW ]i- 
PIV(evanescent wave micro particle image velocimetry) and multilayer nano- 
particle image velocimetry (nPIV)PSl Their accuracy is normally much lower 
than that of force methods due to relatively low optical resolution, system noise 
due to polydispersity of tracers, and difficulties in decoupling flow from diffusion 
(the tracer distribution in the flow field is affected by Taylor dispersion*^). As a 
consequence, it has been always expected that a slippage of the order of a few tens 
nanometers cannot be detected by a velocimetry technique. However, recently 
direct high-precision measurements at the nanoscale have been performed with 
a new optical technique, based on a DF-FCS (double-focus spatial fluorescence 
cross-correlation j^ZEHl (as is schematically shown in Fig. 11.5.1 . As the fluorescence 
tracers are flowing along the channel they are crossing consecutively the two foci, 
producing two time-resolved fluorescence intensities h(t) and hit) recorded in- 
dependently. The time cross-correlation function can be calculated and typically 
exhibits a local maximum. The position of this maximum Tm is characteristic of 
the local velocity of the tracers. Another example of high resolution promising 
applications of FCS consists in determination of average transverse diffusion co- 
efficient to probe slippage.^ Since FCS methods allow consideration of N ~ 10 6 
particles, this gives a satisfactory signal to noise ratio \/JV of order 10 3 , providing 
extremely good resolution as compared with other direct velocimetry methods. 
Coupling with tirfJSO which allows the measurements of the distance of tracers 
from the wall through the exponential decay of an evanescent wave, should fur- 
ther improve the accuracy of approach. 
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Figure 1.5. Schematics of the double-focus spatial fluorescence cross-correlation method. 
Two laser foci are placed along the x axis separated by a distance of a few ftm. They in- 
dependently record the time-resolved fluorescence intensities I\(t) and 12(f)- The forward 
cross-correlation of these two signals yields G(f). Two foci a re sca nned simultaneously 
along the z axis to probe the velocity profile u(z). Adapted froml^TJiS 

In the vicinity of the wall the tracers are submitted to a Taylor dispersion, i.e. 
their diffusion combined with shear enhances the migration speed in the flow di- 
rection. This phenomenon seriously complicates all velocimetry methods since to 
extract the velocity of fluid this effect should be modelled precisely. As shown 
by analysis of DF-FCS dataF 1 large observed values of the apparent slip at the 
hydrophilic wall are normally fully attributed to a Taylor dispersion of nanotrac- 
ers (see Fig. \1.6.i . The data obtained with other velocimetry technique still await 
clarification. We suggest however that some very large values of a hydrophobic 
slip 1 4 ° l 41 l might reflect a Taylor dispersion too. 




Figure 1.6. Typical velocity profile v(z) (open circles) measured in a ~ 100 pirn channel 
with 10~ 4 mol/L NaCl solution. Inset shows the same profile in the vicinity of the wall. 
Error bars were determined from the deviation among repeated measurements. When error 
bars are absent, uncertainties are smaller than the circles. Dots are simulation results. Dash- 
dotted line fits data (the apparent slip length, b a p p ~ 740 mn), and dashed line shows the 
wall location (z = 50.69fim). Adapted frorrP3 
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Many experiments have been performed on the subject, with somewhat con- 
tradicting results. Experimental work focussed mostly on bare (smooth) surface, 
more recent investigations have turned towards rough and structured surfaces, in 
particular super-hydrophobic surfaces. 51 We refer the reader to comprehensive re- 
view article d 6 ! 52 ! for detailed account of early experimental work. In our chapter 
we mention only what we believe is the most relevant recent contribution to the 
subject of flow past 'simple' smooth hydrophobic and rough hydrophilic surfaces, 
which clarified the existing controversies in the field. We focus, however, more on 
the implication of micro- and nanostructuring on fluidic transport, which is still at 
its infancy and remains to be explored. 

1.5. Smooth surfaces: Slippage vs wetting 

From the theoretica l 16 ! 18 ! and simulatiorP^l point of view slippage should not ap- 
pear on a hydrophilic surface, except probably as at very high shear rate.^S A slip 
length of the order of hundred nanometers or smaller is, however, expected for a 
hydrophobic surface! 12 ' 16 ' 20 '^ 

On the experimental side, no consensus was achieved until recently. While 
some experimental data were consistent with the theoretical expectations both for 
hydrophilic and hydrophobic surfaces J 33 l 35 l 56 l some other reports completely es- 
caped from this picture with both quantitative (slippage over hydrophilic surface, 
shear rate dependent slippage, rate threshold for slip, etc) and quantitative (slip 
length of several ptms) discrepancies (for a review see% More recent experiments, 
performed with various new experimental methods, finally concluded that water 
does not slip on s mooth hyd rophilic surfaces, and develops a slip only on hy- 
drophobic surface! 38 ' 47 ' 49 ' 57 ^ One can therefore conclude that a concept of hy- 
drophobic slippage is now widely accepted. 

An important issue is the amplitude of hydrophobic slip. The observed slip 
length reached the range 20-100 nm, which is above predictions of the models 
of molecular slip! 13 ! 19 ! This suggests the apparent slip, such as the 'gas cushion 
model', Eq. lll.2t . Water glides on air, owing to the large viscosity ratio between 
water and air (typically a factor of 50). Experimental values of b suggest that the 
thickness of this 'layer' is below 2 nm. Another scenario of apparent slip such as a 
nanobubble coated surf ac e 22 " 25 ' has to be explored in more details. 

An important conclusion is that it is impossible to benefit of such a nanomet- 
ric slip at separations 0(j.im) and larger, i.e. in microfluidic applications. This is 
why in the discussion of super-hydrophobic slippage below we often ignore a slip 
past hydrophobic solids. However, a hydrophobic slippage is likely of major im- 
portance in nanochannels (highly confined hydrophobic pores, biochannels, etc), 
where ordinary Poiseuille flow is fully suppressed. 

1.6. Rough surfaces 

Only a very few solids are molecularly smooth. Most of them are naturally rough, 
often at a micro- and nanoscale, due to their structure, methods of preparation, var- 
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ious coatings. These surfaces are very often in the Wenzel (impaled) state, where 
solid/liquid interface has the same area as the solid surface (Fig. I1.7."h ). How- 
ever, even for rough hydrophilic Wenzel surfaces the situation was not very clear, 
and opposite experimental conclusions have been made: one is that roughness 
generates extremely large slip,™ and one is that it decreases the degree of slip- 
page.EBEl j± conclusion in favor of slip past rough surfaces is often made based 
on the fact that the hydrodynamic drag force becomes smaller when one of the 
surfaces is rough (see Fig |1.8.t , which is qualitatively consistent with the model 
of slip, Eqs. lll.3l l, (|1.4Jl . Recent experimental data suggested that the description of 
flow near rough surfaces has to be corrected. Such a flow is equivalent to expected 
near a smooth hydrophilic (i.e. no-slip) surface that is located between top and 
bottom of asperities. 38 In other words the correction should be done not for a slip, 
but for a separation from a rough surface. It has also been shown that opposing 
conclusions made byS2HH] simply reflect a different way of a definition of zero sep- 
aration in the AFM (at the top of asperities) and SEA (at their bottom). 




Figure 1.7. Schematic representation of the (a) Wenzel and (b) Cassie pictures with the 
local flow profiles at the gas and solid areas. Reprinted with permission fromPS Copyright 
(2011) by the IOP Publishing Ltd. 



The theoretical description of such a flow represents a very difficult problem. 
Previous theoretical investigations addressed only a case of a far-field flow, where 
spatial variations are small compared to the typical length scale of the problem. 
Most of the articles considered the issue of periodic roughness although some 
recent work exploited non-periodic, spatially homogeneous random roughness.^ 
It has been concluded that these situations can be modeled by the application of 
the Navier condition to an equivalent smooth plane deposited on the top of rough- 
ness |64163 However, there have been also theoretical (far-field) arguments that this 
should be rather a smooth equivalent no-slip wall placed at some distance below 
the top, but still above the bottom of roughness!^ 

This issue was recently resolved in the LB (lattice Boltzmann) simulation 
study^ where the hydrodynamic interaction between a smooth sphere of radius 
R and a randomly rough plane was studied (as shown in Fig. 11.9. \ . Beside its sig- 
nificance as a geometry of SFA/ AFM dynamic force experiments, this allowed one 
to explore both far and near-field flows in a single 'experiment'. The 'measured' 
hydrodynamic force was smaller than predicted for two smooth surfaces (with the 
separation defined at the top of asperities) if the standard no-slip boundary con- 
ditions are used in the calculation. Moreover, at small separations the force was 
even weaker and shows different asymptotics than expected if one invokes slip- 
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Figure 1.8. (a) Hydrodynamic force acting on a rough sphere (R = 20/im) approaching a 
smooth plane. From bottom to top are data (symbols) obtained at the driving speed -6, -10, - 
20, -30, and -40 pm/s. Solid curves show the calculation results obtained for the same speed, 
by assuming no-slip boundary conditions, (b) An AFM image of an apex of a gold-coated 
sphere. Adapted fronP 



page at the smooth fluid-solid interfaces. This can only be explained by the model 
of a no-slip wall, located at an intermediate position (controlled by the density of 
roughness elements) between top and bottom of asperities (illustrated by dashed 
line Fig. ll.7.b ). Calculations based on this model provided an excellent description 
of the simulation data (Fig. |1.9.T >. 

1.7. Super-hydrophobicity and effective hydrodynamic slippage 

On hydrophobic solids, the situation is different from that on hydrophilic solids. 
If the solid is rough enough, we do not expect that the liquid will conform to the 
solid surface, as assumed in the Wenzel or impaled state. Rather air pockets should 
form below the liquid, provided that the energetic cost associated with all the cor- 
responding liquid /vapor interfaces is smaller than the energy gained not to follow 
the solidP This is so-called Cassie or fakir state. Hydrophobic Cassie materials 
generate large contact angles and small hysteresis, ideal conditions for making 
water drops very mobile. It is natural to expect a large effective slip in a Cassie 
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Figure 1.9. Hydrodynamic force acting on a hydrophilic sphere of radius R approach- 
ing a smooth hydrophilic (diamonds), smooth hydrophobic (circles) and randomly rough 
hydrophilic (triangles) wall with <p2 = 4% (adapted frorrp2l). Here F st = 6nt]RU is the 
Stokes drag. The separation h is defined on top of the surface roughness as shown in the 
inset. Simulation results (symbols) compared with theoretical curves: F/F$ t = 1 +9R/ (8h) 
(solid), F/Fst = 1 + 9Rf*/(8h) with /* = f(b/h), calculated with Eq.(TJ) (dash-dotted), 
and F /Fg t = 1 + 9R/ (8[h + s]) (dashed). Values of b and s were determined by fitting the 
simulation data. 

situation. Indeed, taking into account that the variation of the texture height, e, is 
in the typical interval 0.1 — 10 ^m, according to Eq. l|1.2ll we get b = 5 — 500 }im 
at the gas area. The composite nature of the texture requires regions of very low 
slip (or no slip) in direct contact with the liquid, so the effective slip length of the 
surface, b e ff, is smaller than b. Still, one can expect that a rational design of such a 
texture could lead to a large values of b e ff . Below we make these arguments more 
quantitative. 

We examine an idealized super-hydrophobic surface in the Cassie state 
sketched in Fig. I1.7.b where a liquid slab lies on top of the surface roughness. The 
liquid / gas interface is assumed to be flat with no meniscus curvature, so that the 
modeled super-hydrophobic surface appears as a perfectly smooth with a pattern 
of boundary conditions. In the simplified description the latter are taken as no- 
slip (b\ = 0) over solid/liquid areas and as partial slip {bi = b) over gas/liquid 
regions [as we have shown above, b\ is of the orders of tens nm, so that one could 
neglect it since b-i is of the order of tens of }tm]. We denote as 5 a the typical 
length scale of gas/liquid areas. The fraction of solid/liquid areas will be denoted 
(pi = (L — S) I L, and of gas /liquid area </>2 = l — <p\ = 5/L. Overall, the descrip- 
tion of a super-hydrophobic surface we use here is similar to those considered in 
In this idealization, some assumptions may have a possible influence 
on the friction properties and, therefore, a hydrodynamic force. First, by assuming 
flat interface, we have neglected an additional mechanism for a dissipation con- 
nected with the meniscus curvature. 76 78 Second, we ignore a possible transition 
towards impaled (Wenzel) state that can be provoked by additional pressure in the 
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liquid ph ase.M Third, we do not take into account the circulation of air in the 
eas phased™ and some possible other effects such as a depinning of the contact 
lineP 

Finally, for the sake of brevity we focus below only on the canonical microflu- 
idic geometry where the fluid is confined between flat plates, and only on the 
asymmetric case, where one (upper) surface is smooth hydrophilic and another 
(lower) represents a super-hydrophobic wall in the Cassie state. Such a configu- 
ration is relevant for various setups, where the alignment of opposite textures is 
inconvenient or difficult. We also restrict the discussion by a pressure-driven flow 
governed by the Stokes equations: 

rjV 2 u = Vp, V ■ u = 0, (1.5) 

where u is the velocity vector, and p is pressure. Extensions of our analysis to 
study other configuration geometries and types of flow would be straightforward. 

1.7.1. Anisotropic surfaces. 

Many natural and synthetic textures are isotropic. However, it can be interesting 
to design directional structures, such as arrays of parallel grooves or microwrin- 
kles, that consequently generate anisotropic effective slip in the Cassie regime. The 
hydrodynamic slippage is quite different along and perpendicular to the grooves. 
Axial motion is preferred, and such designs are appropriate when liquid must be 
guided. There are examples of such patterns in nature, such as the wings of but- 
terflies or water striders. 




Figure 1.10. Sketch of tensorial hydrodynamic slip. The normal traction f„ exerted by the 
fluid on an anisotropic surface produces an effective slip velocity Au = M f ,, in a different 
direction. At the molecular level, the interfacial mobility tensor M is related to the trajecto- 
ries of diffusing particles in the interfacial region, such as the one shown. Reprinted with 
permission from:?* 5 ' Copyright (2008) by the Cambridge University Press. 

A generalization of the Navier mobility condition is^ 

Au = M ■ (ft ■ tr) (1.6) 

where f„ = n ■ cr is the fluid force (normal traction) on the interface, <r is the local 
stress tensor, and M is an interfacial mobility tensor. As shown in Figure [T. 10.1 the 
effective slip vector is generally misaligned with the force vector for an anisotropic 
surface, which has been analyzed in a number of studies.™ For anisotropic sur- 
faces, the mobility must be a second-rank tensor M = {M,,}. In the case of a 
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Figure 1.11. Sketch of a flat channel of thickness H wall with notation for directions along 
the plates. One wall represents an anisotropic super-hydrophobic texture. Adapted frornQS 



Newtonian fluid, b e ff = M.T], and we get a tensorial version of il.li , as discussed 



(1.7) 



where (u|^) is the effective slip velocity, averaged over the surface pattern and 
n is a unit vector normal to the surface A. The second-rank effective slip tensor 
b e ff = {bjj*} characterizes the surface anisotropy and is represented by symmetric, 
positive definite 2x2 matrix diagonalized by a rotation: 



?eff 




cos v sin V 
— sin 6 cos 8 

J 



(1.8) 



As proven irP^ for all anisotropic surfaces the eigenvalues b^ ff and b^ of the slip- 
length tensor correspond to the fastest (greatest forward slip) and slowest (least 
forward slip) directions, which are always orthogonal (see Fig. ll.lljV 

To illustrate the calculation of the slip-length tensor, below we consider the 
geometry where the liquid is confined between two plates separated by a distance 
H, and one of them represents a super-hydrophobic striped wall (Fig. ll.12Tk ). 

The effective slip lengths in eigendirections (which are in this case obviously 
parallel and orthogonal to stripes) has been recently calculated^ by using the dual 
series technique suggested in work.^lZSI Fig. ll.13lk shows the typical calculation re- 
sults (the numerical example corresponds tob/L = 20 and tp2 = 0.75), and demon- 
strates that the effective slip lengths increase with H and saturate for a thick gap. 
This points to the fact that an effective boundary condition is not a characteristic 
of liquid / solid interface solely, but depends on the flow configuration and inter- 
play between typical length scales, L, H, and b, of the problem. Next we discuss 
asymptotic limits (of small and large gaps) of this semi-analytical solutionP^l 
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Figure 1.12. Special textures arising in the theory: (a) stripes, which attain the Wiener 
bounds of maximal and minimal effective slip, if oriented parallel or perpendicular to the 
pressure gradient, respectively; (b) the Hashin-Shtrikman fractal pattern of nested circles, 
which attains the maximal/minimal slip among all isotropic textures (patched should fill up 
the whole space, but their number is limited here for clarity); and (c) the Schulgasser and 
(d) chessboard textures, whose effective slip follows from the phase-interchange theorem. 
Adapted frorrP 




Figure 1.13. Eigenvalues ir ff (solid curve) and b^ (dashed curve) of the slip length tensor 
for stick-slip stripes of period L with local slip length at liquid-gas interface b/L = 20 and 
slipping area fraction tfo = 0.75 as a function of the thickness of the channel, H. Reprinted 
with permission from. 14 Copyright (2011) by the IOP Publishing Ltd. 

1.7.1.1. Thin channel 

In the limit of H <C L we get 

|| bHcp 2 x bH<p 2 

These expressions are independent on L, but depend on H, and suggest to distin- 
guish between two separate cases. 
If b <C H we obtain 

b^ ~ fejg ~ b<p 2 , (1.10) 

so that despite the surface anisotropy we predict a simple surface averaged effec- 
tive slip. Although this limit is less important for pressure-driven microfluidics, it 
may have relevance for amplifying transport phenomena!^ 
When H <C b we derive 

h l^ u %' ^ii- (i-ii) 
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The above formula implies the effective slip length is generally four times as large 
for parallel versus perpendicular pressure driven flow. Both asymptotic results, 
Eqs. dl.10l ) and dl.lU , are surprising taking into account that for anisotropic Stokes 
flow in a thick channel factor of two is often expected as reminiscent results for 
striped pipesj^ sinusoidal grooves^ and the classical result that a rod sediments 
twice as fast in creeping flow if aligned vertically rather than horizontally.^ A 
very important conclusion from our analysis is that this standard scenario can sig- 
nificantly differ in a thin super-hydrophobic channel, by giving a whole spectrum 
of possibilities, from isotropic to highly anisotropic flow, depending on the ratio 
b/H. 




10~ 2 10~ 1 10° 10 1 10 2 10 3 



(b) b/H 

Figure 1.14. Effective slip length, b e f[/H, versus cp2 [for b/H = 1] (a) and b/H [for 
(p2 = 0.5] (b) in a thin gap limit, H <C L. SH surfaces are: anisotropic stick-slip stripes 
attaining Wiener bounds (dashed curves), isotropic textures attaining Hashin-Shtrickman 
bounds (upper - solid, lower - dash-dotted curves) and satisfying the phase-interchange 
theorem (circles). Reprinted with permission fromiM Copyright (2011) by the IOP Publish- 
ing Ltd. 



Note that in case of a thin channel the flow can be described by an expression 
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of Darcy's law, which relates the depth-averaged fluid velocity to an average pres- 
sure gradient along the plates through the effective permeability of the channel. 
The permeability, K e ff, is in turn expressed through effective slip length b e ff, and 
permeability and slip-length tensors are coaxial. Such an approach allows one to 
use the theory of transport in heterogeneous media, 89 which provides exact results 
for an effective permeability over length scales much larger than the heterogene- 
ity. This theory allows one to derive rigorous bounds on an effective slip length 
for arbitrary textures, given only the area fraction and local (any) slip lengths of 
the low-slip (b\) and high-slip (^2) regions P2I22I These bounds constrain the attain- 
able effective slip and provide theoretical guidance for texture optimization, since 
they are attained only by certain special textures in the theory. In some regimes, 
the bounds are close enough to obviate the need for tedious calculations of flows 
over particular textures. In particular, by using the general result of the theory of 
bounds^ one can easily derive Eq. l|1.9t and its limits, as well as to prove that for a 
thin channel longitudinal (transverse) stripes satisfy upper (lower) Wiener bounds, 
i.e. provide the largest (smallest) possible slip that can be achieved by any texture. 
We remark and stress that according to results 29 stripes should not be necessarily 
periodic. 

Typical dependence of Wiener bounds for b e ff /Hon^ (at fixed b/H) and on 
b/H (at fixed (p 2 ) is shown in Fig. II. 141 which well illustrates that the key parame- 
ters determining effective slip in the thin channel is the area fraction of solid, cf>i, in 
contact with the liquid. If this is very small (or (p 2 — > 1), for all textures the effective 
slip tends to a maximum value, b e ff — > b. We can also conclude that maximizing b 
also plays a important role to achieve large effective slip. 



1.7.1.2. Thick channel 



In the opposite case of infinitely large thickness (H 3> L) the effective slip lengths 

ZQl 
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(1.12) 



(1.13) 



The above results apply for a single surface, and are independent on H. However, 
these expressions for effective slip lengths depend strongly on a texture period 
L. When we again derives the area-averaged slip length, Eq. dl.10l l- When 

b 3> L, expressions l|1.12b and ( 11.13t take form 

2 

that coincides with the result obtained by Lauga and Stone^ for the ideal slip 
(b — > 00) case. We stress that the commonly expected factor of two for the ratio of 
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eff' 



(1.14) 
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bjlg and is predicted only for very large b/L. In all other situation the anisotropy 
of the flow is smaller and even disappears at moderate b/L. 

1 .7.2. Isotropic surfaces. 

As stressed above most solids are isotropic, i.e. without a preferred direction. Un- 
fortunately, from the hydrodynamic point of view, this situation is more compli- 
cated than considered above. Below we discuss only some aspects of the hydro- 
dynamic behavior in the thin and thick channel situations. For a thin channel, 
arguments are based on the already mentioned theory of transport in heteroge- 
neous medial and derived bounds on an effective slip length (effective slip length 
and permeability tensors are now becoming simply proportional to the unit ten- 
sor I) for arbitrary isotropic textures {22I2QI For a thick channel, the only available 
arguments are based on the scaling theory and numerical calculations^ which 
however provide us with some guidance. 

1 .7.2.1 . Thin channel 

If the only knowledge about the two-phase isotropic texture is <p\, <p 2 , then the 
Hashin-Shtrikman (HS) bounds apply for the effective permeability, b y givi ng the 
corresponding upper and lower HS bounds for the effective slip lengthP2E3 These 
bounds can be attained by the special HS fractal pattern sketched in Fig. I1.12.b . 
For one bound, space is filled by disks of all sizes, each containing a circular core 
of one component and a thick ring of the other (with proportions set by the concen- 
tration), and switching the components gives the other bound. Fractal geometry 
is not necessary, however, since periodic honeycomb-like structures can also attain 
the bounds. The general 

solutiorPSQl 

allows to deduce a consequential analytical 
results for an asymmetric case considered in this paper. The upper (HS) bound can 
be then presented as 

bH(p 2 (2H + 5b) 



eff H(2H + 5b) + bfa (5H + b) ' 
and the lower (HS) bound reads 



(1.15) 



_ 2bH<p 2 

h ^-2H + 5b<p{ (L16) 

At small b/H we get Eq. (|1.10t , and at large b/ H these give for upper and lower 
bounds 

8<fr 50i 

correspondingly. 

Finally, phase interchange results^ can be applied for some specific patterns 
(Fig. ll.l2~!b ,d). The phase interchange theorem states that the effective permeability 
K e ff (b\, b 2 ) of the medium, when rotated by 7r/2, is related to the effective perme- 
ability of the medium obtained by interchanging phases 1 and 2, viz. K e g(&2/ b\): 

[R • Kefffa, b 2 ) ■ R f ] ■ K ef£ (b 2 , bi) = K&l 
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where b\ t 2 are the local slip lengths for each medium, R is the rotation tensor and 
R* is its transpose. In the particular case of a medium which is invariant by a re/ 2 
rotation followed by a phase interchange, the classical result follows: 

K e ff = V K 1 K 2 
Obviously <pi — <p2 = 0.5 for such a medium so that: 

b ei( = 3H - H, (1.18) 

! + ■ 



H + b 

At b/H -C 1 we again derive Eq. (|1.10b . indicating that at this limit all textures 
show a kind of universal behavior and the effective slip coincides with the aver- 
age. Thi s sug gests that the effective slip is controlled by the smallest scale of the 
problem ] 16 ! 70 ! so that at this limit b eii is no longer dependent on H, being propor- 
tional to b only. If b/H 3> 1 we simply get 

bes = y (1-19) 

again suggesting a kind of universality, i.e. similarly to anisotropic stripes (cf. 
Eq. ll.lH , in this limit b e # / H for all isotropic textures almost likely scale as cx (p2/<p\- 
The results for these special textures are included in Fig. 11.14.1 which shows 
that Hashin-Strickman bounds are relatively close and confined between Wiener 
ones. For completeness, we give in Table 11.11 a summary of main expressions for 
an effective slip in a thin channel. 

Table 1.1. The effective slip length fc e ff for different textures (shown in Fig. I1.12.| l in a thin gap limit 

(H<L) 

Texture fo e ff 
Wiener upper bound (longitudinal stripes) bH<p2 



H + bep! 

Wiener lower bound (transverse stripes) — bHfo — 
Hashin-Shtrickman upper bound (Hashin-Shtrickman fractal, honeycomb-like texture) 



H(2H + 5b) + fo^ 1 (5H + fc) 

Hashin-Shtrickman lower bound (Hashin-Shtrickman fractal, honeycomb-like texture) 2H~~5^ — 

3M 

Phase interchange patterns (Schulgasser texture, family of chessboards) — H 

4- Jl + 3b/(H+b) 



1.7.2.2. Thick channel 

For this situation the exact solution was not found so far. Nevertheless, some 
simple scaling expressions have been proposed for a geometry of pillarsj^ES] D y 
predicting b e {{ c< L/(7iy / ^7) (cf. scaling results for stripes b e {{ c< L / / (p\)) . 
This simple result would deserve some analytical justification, which has not been 
performed up to now, despite some recent approximate analysis!^ Recent semi- 
analytical and numerical results confirm this scaling dependence,^ which, in par- 
ticular, suggests that in a thick channel the array of pillars will give larger effective 
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slip than longitudinal stripes for a sufficiently small solid fraction. This is opposite 
to a prediction for a thin channel, where an array of longitudinal stripes provides 
the largest possible slip that can be achieved by any texture, whether isotropic or 
anisotropic. 



1.8. Other special properties of super-hydrophobic surfaces 

As we see above, hydrophobic Cassie materials generate large and anisotropic ef- 
fective slippage compared to simple, smooth channels, ideal situation for various 
potential applications. A straightforward implication of super-hydrophobic slip 
would be the great reduction of the viscous drag of thin microchannels (enhanced 
forward flow), and some useful examples can be found inP Below we illustrate 
the potential of super-hydrophobic surfaces and possibilities of the effective slip 
approach by discussing a couple of other applications. Namely, we show that 
optimized super-hydrophobic textures may be successfully used in a passive mi- 
crofluidic mixing and for a reduction of a hydrodynamic drag force. 



1 .8.1 . Transverse flow 



The effective hydrodynamic slip 1 3 ! 28 ! 30 ! of anisotropic textured surfaces is generally 
tensorial, which is due to secondary flows transverse to the direction of the applied 
pressure gradient. In the case of grooved no-slip surfaces (Wenzel state), such 
a flow has been analyzed for small height variations^ and thick channels,^ and 
herringbone patterns have been designed to achieve passive chaotic mixing during 
pressure-driven flow through a microchannell^H 

It has recently been demonstrated that similar effects may be generated by a 
super-hydrophobic Cassie surface. The transverse flow due to surface anisotropy 
is generated only in the vicinity of the wall and disappears far from it, which has 
already been observed in experiment.^ In other words, the effective velocity pro- 
file is 'twisted' close to the super-hydrophobic wall (see Fig. ll.l5.k ). 

To evaluate the transverse flow the velocity profile has been integrated across 
the channel to obtain the flow vector: 

H 

Q=/(u(z))dz (1.20) 



with the components: 

o-H 3 
7 12" 
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(1.21) 



Qv = - 



n 4 



sin u cos ( 



(H + 4)(H- 

To optimize the texture, channel thickness and the angle 8 between the directions 
of stripes and the pressure gradient, so that \Qy/ Qx\ is maximum providing the 
best transverse flow. 



^ff) 



(1.22) 
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Figure 1.15. (a) Scheme of a transverse flow generation, (b) Fraction of flow vector Q 
components (maximized over 8) as function of channel thickness for local slip b/L = 1000 
and slip fraction (fo = 0.5 (solid line), 0.2 (dashed) and 0.9 (dash-dot). Reprinted with 
permission from; 14 Copyright (2011) by the IOP Publishing Ltd. 



The maximization in respect to 9 indicates that the optimal angle isj^ 



± arctan 



l+4fr; ff /H)(l + fr e yH) 

L(i+z4/H)(i+4b e yji)j 

and the value of the maximum reads 
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2 V tan# max 



(1.23) 



(1.24) 



We conclude, therefore, that since H is fixed, the maximal \ Qy/ Q x \ corresponds to 
the largest physically possible b, i.e. the perfect slip at the gas sectors. 

To optimize the fraction of the slipping area, <p2, we should now exploit results 
for effective slip lengths obtained above. Fig U.15.b shows the computed value 
°f iQy/Qxl vs - H/L for several cp2- The calculations are made using the value 
of 9 defined by Eq. dl.23l l. The data suggest that the effect of <p2 on a transverse 
flow depends on the thickness of the channel. For a thick gap the increase in gas 
fraction, <p2, augments a transverse flow. This result has a simple explanation. For 
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an infinite channel b^'J~/H <C 1 (see Fig. (|1.14.b )), which gives 



Qv 



3Ah * (1.25) 



i.e. in a thick channel the amplitude of a transverse flow is controlled by the dif- 
ference between eigenvalues of the effective slip tensor, Ab e y = b^ ff — b^, which 
increases with <p2 as follows from the above analysis. We stress however, that since 
\Qy/ Qx\ H~^, the mixing in a thick super-hydrophobic channel would be not 
very efficient. A more appropriate situation corresponds to a thin channel as it is 
well illustrated in Fig ll.l5.b . We see, that the largest transverse flow can be gen- 
erated at intermediate values of (p2- The limit of thin channel has recently been 
studied in details by using a general theory of mathematical bounds]^ and the 
optimum value of <p2 = 0.5 corresponding a numerical example in Fig ll.l5.b has 
been rigorously derived.^ 

An important conclusion from the above analysis is that the surface textures 
which optimize transverse flow can significantly differ from those optimizing ef- 
fective (forward) slip. It is well known, and we additionally demonstrated above, 
that the effective slip of a super-hydrophobic surface is maximized by reducing 
the solid-liquid area fraction <p\ . In contrast, we have shown that transverse flow 
in super-hydrophobic channels is maximized by stripes with a rather large solid 
fraction, (pi = 0.5, where the effective slip is relatively small. 



1 .8.2. Hydrodynamic interactions 




x 



Figure 1.16. Sketch of a hydrophilic disk approaching super-hydrophobic striped disk. 

As a consequence of the superlubrication potential, a hydrophobic texture 
could significantly modify squeeze film drainage between surfaces. It is of ob- 
vious practical interest to minimize the hydrodynamic resistance, F, to approach 
of surfaces. 

For a Reynolds problem, where a disk of radius R is driven towards (in our 
case) a super-hydrophobic plane with a velocity U (see Fig. ll.16l we should solve 
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Eqs.l ll.5t by applying the effective tensorial boundary condition, Eq. (L7) , at the 
super-hydrophobic wall. This allows to derive a general expression for hydrody- 
namic force satisfying the condition p = po at the edge of the disll^ 

3 tztjUR 4 

F = 2 H 3 fell = F R/eff' I 1 - 26 ) 

where Fr represents the classical solution of creeping flow equations of the 
Reynolds lubrication theory,® and the correction for an effective slip is 



F 

/eff — ~F~ 
* R 



H + 4jg(H) H + 4b^(Hy ~ 
H + bl(H) H + ^(H)_ 



(1.27) 



-effv 

Thus the effective correction for a super-hydrophobic slip is the harmonic mean of 
corrections expressed through effective slip lengths in two principal directions, 

/e*f f = 2 (4j + ^i) a- 28 ) 

V/eff ■'eff / 

In case of isotropic textures, all directions are equivalent with b^ ff = = b e ff, 
so we get 

r F _ H + MH) 

Fr ~ H + 4b e ff(H) (L29) 

Obviously, the case b^ = b^- fi = corresponds to /* ff = 1 and gives the Reynolds 
formula. 

The expression for /* ff is very general and relates it to the effective slip length 
of the super-hydrophobic wall and the gap. In order to quantify the reduction of 
a drag force due to a presence of a super-hydrophobic wall, this expression can 
be used for all anisotropic and isotropic textures, where analytical or numerical 
expressions for b^ have been obtained. 

An important consequence of Eq. l|1.27|l is that to reduce a drag force we need 
to maximize the ratio b e n / H, but not the absolute values of effective slip itself. 
This is illustrated in Fig. I1.17.a . where values presented in Fig. I1.13.b were used to 
compute the correction for effective slip, /* ff as a function of the gap. At small H/L 
our calculations reproduce the asymptotic values predicted by Eqs. (1.11) . They 
however vanish at large distances, where b^~ / H are getting negligibly small. The 
useful analytical expressions for f* {{ corresponding to a configuration of stripes are 
presented in Table [L2l Similar estimates for the most important situation of a thin 
gap can similarly be done for some isotropic textures, and we include these results 
into Table[L3] 

The results presented in Tables 11.21 and 11.31 show that the key parameter de- 
termining reduction of drag is the area fraction of gas, <p2, in contact with the liq- 
uid. This is illustrated in Fig. I1.17.b , where (using a relatively large b/H) Hashin- 
Strickman bounds for /* {f are plotted versus (p2- If this is very small (or (p\ — > 1) for 
all textures, the correction for slip tends to its absolute maximum, /* ff — » 1 . In the 
most interesting limit, (p2 — > 1, we can achieve the minimum possible value of cor- 
rection for effective slip, /* ff — »■ 1 / 4 provided b/H is large enough. We also stress 
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that the results for stripes are confined between Hashin-Strickman bounds for . 
In other words, isotropic textures might be the best candidates for a reduction of a 
drag force. 

Another important point to note would be that there exists no universal re- 
lationship between and the effective Cassie angle, cosflj , of the super- 
hydrophobic surface. Indeed, according to the Cassie equation, cos 8^. = (1 — 
(j>2) cos 6 — (p2^ where cos 6 is the thermodynamic contact angle at the smooth sur- 
face. It is well seen already from Fig. I1.17.b that the dependence of /* ff on cp2 is 
much more complex. In particular, at fixed H it depends on texture. In general, 
this also reflects the fact that there is no universal relationship between b eff and <p2, 
despite some recent attempts to establish itf^l 

Table 1.2. Asymptotic expansions for 
the force correction factor f*„ in case of 
a striped surface. 
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4 32 L ln(sec(7r^ 2 /2)) 
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H«min{L, b} ^<h) 
8 + 9^2-9^| 



Table 1.3. Correction factor /* ff for some specific isotropic 
patterns in a thin gap limit. 

Pattern li«H«I H < min{L, b} 

Hashin-Strickman upper bound 1 — ^ 



Hashin-Strickman lower bound 
Phase interchange textures 1 



H 5 + 302 

_ 3fc02 8-3^2 



H 4(2 + 30) 

3b_ 1 

2H 2 



1.9. Interfacial transport phenomena 
1 .9.1 . Hydrophobic surface 

Besides this drag-reduction potential in pressure-driven flows, it has recently been 
predicted that hydrophobic slippage could be best exploited in surface-driven 
transport. 

Electro-osmosis (EO), i.e. flow generation by an electric field, originally devel- 
oped in colloid science, is currently experiencing a renaissance in microfluidics. 
It may be considerably amplified by surface slippage, even for nanometric slip 
length. The reason for this amplification is that the electric double layer (EDL), 
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(a) 




0.2 0.4 0.6 0.8 1 



(b) ^ 

Figure 1.17. (a) The correction factor /* ff = F/Fr for hydrodynamic resistance force ex- 
erted on disk interacting with super-hydrophobic stripes vs. dimensionless gap width H/ L. 
Solid curves correspond to local slip length b/L = 10 (from top to bottom <p2 = 0.2, 0.5 and 
0.9), dashed curves - to b/L = 0.1 (from top to bottom (p2 = 0.2 and 0.5), dash-dotted curve 
- to b/L = 0.01 and <p2 = 0.5. (b) The plot of /* ff versus <p2 for a thin gap (H <C L) and sev- 
eral super-hydrophobic patterns: anisotropic stripes (solid line), isotropic textures attaining 
Hashin-Strickman bounds (dashed and dash-dotted lines) and isotropic Schulgasser struc- 
ture (diamond), all with local slip b/H = 10. Reprinted with permission fromU^ Copyright 
(2011) by the IOP Publishing Ltd. 

characterized by the Debye screening length = k _1 defines an additional length 
scale of the problem comparable to b. According to classical formula!- 100 ! 101 ! 

u = -^(l + k) = -^(l + k) (1.30) 

where u is the electro-osmotic velocity (outside of the double layer), qo is the sur- 
face charge density, £ = qo/ tee is the zeta-potential across the diffuse (flowing) 
part of the double layer, e is the permittivity of the solution, and Eq is the tan- 
gential electric field. Therefore, the flow can potentially be enhanced for a thin 
compared to b EDL, i.e. when k!) > 1 (see Fig. 11.18. L which has been proven 
experimentallyP^HSU 
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(b) 



Figure 1.18. Sketch of the influence of hydrophobic (a) and superhydrophobic (b) slippage 
on the electro-osmotic flow. 



This expression can be generalized to other interfacial transport phenomena, 
such as diffusio-osmosis and thermo-osmosis, as it has been well discussed in! 16 l 86 l 
These correspond to the induction of a flow by the gradient of a solute concentra- 
tion for the former and by a gradient of temperature for the latter!i22l The flow ve- 
locity for these two important phenomena is proportional to the applied gradient 
of concentration or temperature. Similarly to electro-osmosis, for both, slippage 
amplifies the velocity with the factor (1 + b/X), where A is a thickness of a thin 
interface layer (of the order of the range of interaction of the solute with the solid 
surface). 



1 .9.2. Super-hydrophobic surface 

It is now attractive to consider electro-osmotic flow over super-hydrophobic sur- 
faces, whose texture can significantly amplify hydrodynamic slip. Eq. JT.30b sug- 
gests that a massive amplification of EO flow can be potentially achieved over 
super-hydrophobic surfaces. However, the controlled generation of such flows is 
by no means obvious since both the slip length and a charge distribution are in- 
homogeneous and anisotropic. Despite its fundamental and practical significance, 
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electro-osmotic flow over super-hydrophobic surface has so far received little at- 
tention. Only recently such a flow has been investigated past inhomogeneously 
charged slipping surface in the case of a thick channel (H 3> L) and perfect slip 
(b — s- oo) at the gas sectors.ESllIQil 




Figure 1.19. The electro-osmotic slip velocities (thin EDL) for eigendirections of super- 
hydrophobic stripes. Dash-dotted line corresponds to the case of uncharged gas sectors 
<?1 — <?0> <?2 = 0: transverse and longitudinal EO flows coincide with each other. Upper 
solid curve corresponds to longitudinal, and upper dashed - to transverse EO flow for a 
surface of a constant charge q\ = qi = <?o- Lower solid curve corresponds to longitudinal, 
and upper dashed - to transverse EO flow for a surface with c\\ = —c\i = —qg. 

The general result for a thin EDL (jcL 3> 1) has been formulated irP^Sl 

En 

u = (qil + q 2 b ei {K) (1.31) 

t]K 

by using the Lorentz reciprocal theorem for the Stokes flow. Here I is unity tensor, 
Cj\ and c\2 are me surface charge density at the solid and gas regions, correspond- 
ingly. This expression shows that surface anisotropy generally leads to a tenso- 
rial electro-osmotic mobility To highlight the effect of anisotropy, we focus now 
on patterned super-hydrophobic surface consisting of periodic stripes, where the 
surface charge density varies only in one direction, likewise the slip lengths. By 
combining Eq. il.31t with Eqs. lll.14t we get general results and important limiting 
situations for stripes summarized in Table [T~4l and discussed below. 

In case of uncharged gas interface (q\ = cjo, q 2 — 0) we predict the simple 
Smoluchowski formula. 103 In other words, there is no EO flow enhancement, and 
the flow is isotropic despite anisotropy of a surface (see Fig. II .191 . This surprising 
result has been confirmed by molecular dynamic simulations^^ and later analy- 
sis IIOJJ ]\j (- e ma t f or uncharged gas interface with thick EDL (kL <C 1), the results 
are qualitatively different, and for uncharged gas interface the EO flow remains 
tensorialP^J However, since kL is small and since the electro-osmotic velocity is 
proportional to the fraction of a charged area,<p\, electro-osmotic flow becomes 
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Table 1.4. Electro-osmotic slip past super-hydrophobic stripes in a thin EDL limit 



'cff 



General case 
<7i = <?o< <?2 = 
9i = <?2 = ?o 
<?i = -<?2 = 90 



E 7 kL 
— <7i +</2 — I" 




7T(* 2 



2 

2 



?l +<?2 — In 

t]K 1 




suppressed (compared to predicted by the Smoluchowski equation) despite a large 
effective slipP^ 

In case of a charged gas interface (only) a considerable enhancement of electro- 
osmotic flow is possible. For a uniformly charged (qi = qi = qr>) anisotropic super- 
hydrophobic surface the expression for electro-osmotic flow can be transformed to 



r]K 



(I + b eff K) , 



(1.32) 



which might be seen as a tensorial analog of Eq. ^1.30t . Fig. 11.191 includes theo- 
retical results calculated with Eq. lll.32|l for a geometry of stripes, and is intended 
to demonstrate that the flow is truly anisotropic and can exhibit a large enhance- 
ment from effective hydrodynamic slip, possibly by an order of magnitudes. We 
stress that such an enhancement is possible even at a relatively low gas fraction, i.e. 
when b e ff is relatively small (but the amplification ratio, b e ff /Ad, might be huge). 
Obviously, in case of isotropic super-hydrophobic surface, Eq. l ll.32|l transforms to 



t]K 



effK) , 



(1.33) 



and the amplification of EO flow at the isotropic super-hydrophobic surface (such 
as observed in a recent experimental) might serve as a very strong evidence in 
favor of a charge at the liquid-gas interface. 

An interesting scenario is expected for oppositely charged solid and gas sectors. 
In this case Eq. l ll.31|l transforms to 



r]K 



(I - beffK) . 



(1.34) 



The calculation results for this situation are also included into Fig. 11.19.1 and sug- 
gest a very rich fluid behavior. We see, in particular, that inhomogeneous surface 
charge can induce EO flow along and opposite to the field, depending on the frac- 
tion of the slipping area. Already a very small fraction of the gas sectors would 
be enough to reverse the effective EO flow. Another striking result is that electro- 
neutral surface ((q) = (p\q\ + fatfl = 0) can generate extremely large EO slip. With 
our numerical example this corresponds to <p2 = 0.5. 

Of course, the area of research connected with interface transport phenomena 
is still at its infancy. Thus, the role of the surface conductance just started to be 
probedPa 

Beside that, many assumptions exploited above should obviously be 
relaxed. For example, in future, we suggest as a fruitful direction to consider 
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electro-osmotic flow in a thin gap and by assuming a partial slip at the gas sec- 
tors. It will be very important to investigate transverse electrokinetic phenomena, 
that could be greatly amplified by using striped super-hydrophobic surfaces. Note, 
that if the charge is varied along the direction of electric field, the fluid close to a 
super-hydrophobic wall is pulled periodically in opposite directions. As a result, 
the recirculation rolls should develop on a scale proportional to a texture size, L. 
This should provide an additional opportunity for a mixing, similar to described 
inPlbut hopefully much faster. 

1.10. Conclusion and future directions 

With recent progress in micro- and nanofluidics new interest has arisen in deter- 
mining forms of hydrodynamic boundary conditions! 3 ! 28 ! 30 ! In particular, advances 
in lithography to pattern substrates have raised several questions in the modeling 
of the liquid motions over these surfaces and led to the concept of the effective 
tensorial slip. These effective conditions capture complicated effects of surface 
anisotropy and can be used to quantify the flow over complex textures without 
the tedium of enforcing real inhomogeneous boundary conditions. 

This chapter has discussed the issue of boundary conditions at smooth hy- 
drophobic and rough hydrophilic surfaces, and has then given the especial empha- 
sis to the effective boundary conditions for a flow past hydrophobic solid surfaces 
with special textures that can exhibit greatly enhanced ('super') properties, com- 
pared to analogous flat or slightly disordered surfaces. Research on these super- 
hydrophobic materials during past decades has mostly focused on their extreme 
non-wettabilityPSl However, now the field has moved beyond wetting towards 
transport phenomena! 14 ! 109 ! An effective slip becomes the main parameter to quan- 
tify the effective transport properties near super-hydrophobic surfaces, which is in 
contrast to a traditional approach, based on the use of the effective contact angle. 

We have discussed formulas describing effective boundary conditions for 
pressure-driven flow past super-hydrophobic textures of special interest (such as 
stripes, fractal patterns of nested circles, chessboards, and more). The predicted 
large effective slip of super-hydrophobic surfaces compared to simple, smooth 
channels can greatly lower the viscous drag of thin microchannels. 

The power of the tensor formalism and the concept of effective slippage has 
then been demonstrated by exact solutions for two other potential applications: 
optimization of the transverse flow and analytical results for the hydrodynamic 
resistance to approach of two surfaces. These examples demonstrate that prop- 
erly designed super-hydrophobic surfaces could generate a very strong transverse 
flow and significantly reduce the so-called 'viscous adhesion'. Finally, we have 
discussed how super-hydrophobic surfaces could amplify electrokinetic pumping 
in microfluidic devices. 

As we have shown, a combination of wetting and roughness provides many 
new and very special hydrodynamic properties of surfaces, which could be ex- 
plored with the formalism discussed here. This should allow the local slip ten- 
sors to be determined by global measurements, such as the permeability of a tex- 
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tured channel as a function of the surface orientations or the hydrodynamic force 
exserted on the body approaching a super-hydrophobic plates. They may also 
guide the design of super-hydrophobic surfaces for microfluidic lab-on-a-chip and 
other applications. 

Despite the impressive advances in hydrodynamic slippage phenomena, many 
challenges remain, both theoretical and experimental, fundamental and practi- 
cal. In particular, there is quantitative discrepancy between theory and experi- 
ment (measured effective slip are persistently lower than the theoretical predic- 
tions). This discrepancy reflects gaps in our fundamental understanding of flow 
past rough and hydrophobic surfaces. Below we briefly discuss promising future 
directions for research, both experimental and theoretical. 

• Experiments in thick 

channelsSMHH 

have established that hydrodynamic 
flows are generally slower than one would expect from theoryP^ Current analyti- 
cal models of super-hydrophobic effective slip are based on the idealized model of 
a heterogeneous surface with patches of boundary conditions, and mostly neglect 
a number of dissipation mechanisms in the gas phase and at the interface. The 
effects associated with different aspects of the gas flow and meniscus curvature 
must be included in the models. Regardless recent semi-analytical and numerical 
analysi o 78 ! 81 ! the goal should remain to find simple analytical formulas, with as few 
adjustable parameters as possible, to fit experimental data. 

• Drainage experiments conducted in the AFM and SFA, which has the abil- 
ity to probe fluid films of nano- and molecular thicknesses, may be able to yield 
detailed information on super-hydrophobic slip. Some data are already avail- 
abiejllMIll] ]-, u t we are unaware of any previous work that has studied how the 
squeeze film drainage between curved surfaces would be modified by the occur- 
rence of the effective slip. Recent theory of hydrodynamic interaction between 
disks^l shed some light on what could happen qualitatively, but cannot be used 
for a quantitative analysis of the hydrodynamic data obtained with the AFM and 
SFA. The same remark concerns the use of a theory of a film drainage between 
smooth hydrophobic surfaces, 12 which is not fully applicable to quantify a super- 
hydrophobic slip. We believe that a challenge for a theory would be to develop a 
theoretical modeling of experimentally relevant sphere vs. plane geometry. This 
will open many possibilities for new experiments, and could revolutionize the 
field. 

• There are many opportunities to design new experiments and to develop 
improved theoretical models for electro-osmotic flow past rough, hydrophobic and 
super-hydrophobic surfaces! 103 ! 104 1 For example, a systematic study of the effect of 
surface texture, amplitudes of local slip lengths, and a role of a channel thickness 
would be interesting. It would likewise be interesting to work with conducting, 
but hydrophobic or super-hydrophobic surfaces. The study of super-hydrophobic 
Cassie surfaces will naturally pose fundamental questions about a mechanism of 
electro-osmosis at the charged gas interface, where both adsorbed ions and their 
screening clouds are mobileP^H 

• Very promising directions are certainly the diffusio-osmosis and thermo- 
osmosis, where solvent flows are induced by gradients of solute concentration or 
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temperature. Being combined with hydrophobic^ or super-hydrophobic slippage, 
these effects could lead to a giant amplification of flows in microchannels, even 
if the liquid-gas interface is uncharged as shown in recent simulation studyPSH 
The combination of these two strategies, i.e. diffusio-/ thermo-osmosis and super- 
hydrophobicity has to be studied theoretically and experimentally, and we expect 
a significant expansion into this very interesting area of research. 
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